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ABSTRACT* 


A  method  of  solution  to  the  buckling  problem  for  shallow 
viscoelastic  shells  is  presented.  The  arches  are  considered 
to  be  pin-ended  at  rigid  supports,  of  arbitrary  initial  shape 
but  of  uniform  cross  section,  and  made  of  an  arbitrary  linear 
viscoelastic  material.  The  lateral  loading  may  be  an  arbitrary 
function  of  both  position  and  time. 

Specific  solutions  are  given  for  several  arches  to  point 
out  the  influence  of  the  material  and  of  asymmetries  in  the 
initial  arch  shape  on  the  deformation  and  buckling  time  of 
the  system.  Two  types  of  viscoelastic  material,  the  Maxwell 
Fluid  and  the  three  parameter  solid,  are  included  in  the  exam¬ 
ples. 

The  viscoelastic  property  of  the  material  in  the  analysis 
is  expressed  in  integral  operator  form,  which  is  shown  to 
lead  to  simple  numerical  time  integrations  in  determining  the 
arch  deflection,  as  well  as  maintaining  the  form  of  the  gov¬ 
erning  equations  that  determine  stability  in  the  same  general 
form  as  the  equations  of  an  elastic  arch.  The  method  of  so¬ 
lution  takes  into  account  the  fact  that  the  arch  may  buckle 
in  either  a  symmetrical  or  an  unsymmetrical  mode. 

It  is  seen  that  arches  composed  of  a  viscoelastic  solid 
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Research,  V/ashington,  D.  C.  Reproduction  in  whole  or  in  part 
is  permitted  for  any  purposes  of  the  United  States  Government. 
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material  may,  if  the  load  is  not  too  great  and  constant  for 
time  large,  reach  a  final  unbuckled  position  as  the  time 
becomes  very  large.  If  this  is  the  case  then  it  is  shown 
that  the  time  integration  that  previously  required  numerical 
integration  can  now  be  approximated  analytically.  Several 
examples  are  given. 
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Chapter  I  Introduction  and  Summary 


Whenever  a  structure  or  structural  member  is  subject 
to  compressive  loads,  the  possibility  of  an  instability  type 
of  failure  must  be  considered.  Many  structures,  when  loaded 
to  a  critical  state,  will  undergo  a  marked  change  in  the 
magnitude  and  character  of  their  deformation,  which  is  not 
the  result  of  any  failure  or  alteration  in  the  mechanical 
properties  of  the  material.  Such  a  change  generally  occurs 
because  one  mode  of  deformation  becomes  unstable  and  the 
structure  deflects  to  another,  stable  mode.  This  change  is 
termed  buckling,  and  the  loading  associated  with  this  criti¬ 
cal  state  is  most  generally  called  the  buckling  load.  How¬ 
ever,  there  are  many  structures  for  which  a  well  defined 
buckling  load  does  not  exist.  For  these  structures  the  use¬ 
ful  load  carrying  ability  must  be  based  on  other  factors, 
such  as  maximum  allowable  stresses  or  deflections. 

Structures  that  are  made  of  materials  which  do  not  ex¬ 
hibit  any  time  effects  have  a  buckling  load  that  is  inde¬ 
pendent  of  time  provided  the  load  is  applied  slowly  enough 
so  that  the  response  can  be  considered  as  static.  However, 
viscoelastic  structures,  because  the  material  exhibits  time 
or  memory  effects  in  that  the  strain  is  a  function  of  the 
stress  history,  may  have  a  buckling  load  that  is  a  function 
of  time  even  for  a  quasi-static  analysis  of  the  problem. 
Moreover,  for  structures  that  exhibit  no  well  defined  buck¬ 
ling  action,  viscoelastic  material  properties  can  greatly 
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influence  their  load  carrying  ability,  especially  if  the  de¬ 
formations  tend  to  reduce  the  stiffness  of  the  structure. 

Buckling  occurs  most  often  in  structures  that  have  at 
least  one  dimension  small  in  comparison  to  the  other  two; 
for  massive  structures  not  falling  into  this  class,  material 
failure  more  commonly  provides  the  critical  condition  on  the 
magnitude  of  loading.  An  explanation  of  why  this  is  so  can 
be  seen  by  investigating  the  relationship  between  elastic 
stability  theory  and  linear  elasticity  theory.  In  elasticity 
theory,  the  equations  of  equilibrium  are  written  for  the  body 
in  its  undeformed  state.  If  the  body  is  massive  it  cannot 
undergo  large  deformations  without  large  strains,  and  since 
elastic  strains  for  most  engineering  materials  are  very  small 
it  is  concluded  that  the  deformations  must  be  very  small  and 
can  be  neglected  in  writing  the  equations  of  equilibrium. 

This  approach  then  leads  to  the  linear  equations  of  elasticity 
which  have  a  unique  solution,  and  as  such  cannot  describe  a 
stability  problem.  By  contrast,  in  stability  theory  the  e- 
quations  of  equilibrium  must  be  written  for  the  body  in  its 
deformed  state.  Blender  bodies  can  undergo  large  deformations 
without  producing  large  strains,  and  consequently  the  defor¬ 
mations  must  be  accounted  for  when  v/riting  the  equilibrium 
equations,  since  they  are  no  longer  restricted  to  being  small. 

It  is  possible  to  distinguish  several  different  kinds 
of  buckling,  depending  on  the  type  of  structure  and  the  man¬ 
ner  of  loading.  By  far  the  most  widely  studied  type  is  the 
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so-called  classical  buckling,  an  example  of  which  is  the 
centrally  loaded  straight  column.  A  second  type  is  snap- 
through  buckling,  which  occurs  in  structures  whose  resistance 
to  an  increase  in  load  decreases  as  the  loading  or  the  de¬ 
formations  increase.  An  example  of  this  is  the  laterally 
loaded  shallow  arch,  the  analysis  of  which,  for  the  visco¬ 
elastic  case,  forms  the  body  of  this  dissertation. 

The  classical  buckling  problem  of  the  centrally  loaded 
straight  column,  as  well  as  the  problem  of  the  initially  im¬ 
perfect  column  which  has  no  well  defined  buckling  load,  has 
been  studied  in  detail  for  both  elastic  and  viscoelastic 
materials  (see  references  [1]  through  [9])-  Because  there 
are  essential  differences  in  the  responses  between  an  elastic 
and  a  viscoelastic  column,  the  results  of  some  of  the  above 
references  will  be  outlined  here. 

Consider  a  centrally  loaded  straight  elastic  column  as 
shown  in  Pig.  1.1a.  As  the  load  P  increases  the  column 
remains  perfectly  straight  until  P  =  P  .  1  ?  P  =  P  , 

cr  Cl* 

the  column  may  deflect  laterally  as  shown  in  Pig.  1.1b,  and 
any  value  of  6  is  an  equilibrium  position.  P^„  is  then 
the  critical  load  at  which  a  change  in  deformation  may  occur. 


^Numbers  in  square  brackets  (]  refer  to  the  list  of 
References . 
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Pig .  1.1 

The  solid  curves,  OAB  and  OAC  of  Fig.  1.1c  show  the 
load  vs .  lateral  deflection  parameter  6  ,  where  6  is  the 
value  of  the  deflection  at  the  midheight  of  the  column.  At 
the  point  A  there  exist  two  different  equilibrium  paths  of 
deformation.  This  condition  is  referred  to  as  a  bifurcation 
of  the  equilibrium  deformations,  and  is  characteristic  of  the 
classical  type  of  buckling,  although  other  kinds  of  buckling 
may  exhibit  the  same  phenomenon.  is  defined  as  the 

classical  buckling  load  for  this  structure,  and  is  given  by 
[1]  as 


2 

tt^EI 


(1.1) 


where  E  is  Young's  modulus  of  the  elastic  material,  I  the 
moment  of  inertia  of  the  cross  section  (assumed  independent 
of  X  ),  and  L  the  column  length. 

_4_ 


If  it  is  assumed  that  the  column  has  some  initial  im¬ 
perfection,  v/hich  is  a  much  more  realistic  approach  than 
assuming  a  perfect  column,  then  the  load  vs.  lateral  deflec¬ 
tion  parameter  6  is  given  by  the  dashed  curve  0*A‘C’  in 
Pig.  1.1c.  Let  yQ(x)  be  a  measure  of  the  initial  imper¬ 
fection,  then 

y^(x)  =  y(x)  when  P  =  0  .  (1.2) 

As  y^Cx)  0  the  load  deformation  curve  0*A'C*  approaches 
the  curve  OAC  but  never  OAB  ;  thus  for  a  non-perfect 
column  the  load  deformation  relation  always  has  a  unique 
one-to-one  correspondence.  There  is  no  bifurcation  point 
and  no  well  defined  buckling  load.  However,  as  can  be  seen 
from  Pig.  1.1c,  the  classical  buckling  load  P  of  the 

w  X 

perfect  column  serves  a  very  useful  purpose  for  the  non¬ 
perfect  column,  in  that  as  P  P^^  the  ratio  6/6^  in¬ 
creases  very  rapidly.  Thus  P  may  be  considered  as  an 
upper  limit  of  useful  loads  that  can  be  supporter  by  the  non¬ 
perfect  column.  As  will  be  shown,  the  same  cannot  be  said 
about  the  critical  buckling  load  of  a  perfect  viscoelastic 
column,  whose  analysis  follows. 

Consider  a  centrally  loaded  straight  viscoelastic  column 
loaded  as  shovm  in  Pig.  1.1a.  If  the  column  does  not  buckle, 
the  only  deformation  is  a  shortening  of  the  column  due  to  the 
uniform  axial  load.  If  buckling  is  to  take  place,  the  lateral 
motion  must  be  governed  by  the  instantaneous  modulus  of  the 
material  since  there  is  no  motion  in  this  direction  before 
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bO^kling.  Based  on  this  the  buckling  load  for  a  perfect 
viscoelastic  column  is 

P-„  =  (1.3) 

cr 

where  E(0)  is  the  instantaneous  elastic  modulus  of  the 
viscoelastic  material,  and  is  defined  as  the  modulus  govern¬ 
ing  the  stress  immediately  after  the  material  has  been  sub¬ 
jected  to  a  suddenly  applied  strain. 

If  now  a  non-perfect  viscoelastic  column  is  considered, 
the  lateral  deflection  is  given  as  a  function  of  time,  ini¬ 
tial  imperfection,  and  load.  The  general  features  of  the 
response  of  a  non-perfect  column  to  a  load  that  is  suddenly 
applied  and  then  held  constant  can  be  brought  out  by  a 
quasi-static  linear  analysis  of  a  column  whose  initial  im¬ 
perfection  consists  of  a  lateral  displacement  in  the  form  of 
a  half  sine  wave. 

For  the  column  as  chovm  in  Fig.  1.2a,  let 


sin 


TTX 


O'”  '  'O  ”  L 
be  the  displacement  before  the  load  is  applied,  and 


(1.4) 


•rr/v  4-  ^ 

Jf  V  / 


x  /  \  —  j  ^  TTX. 

/  s>  Xll 


/  “»  r-  \ 


the  total  displacement  after  loading. 
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Prom  the  solution  given  in  [2]  for  several  different 
viscoelastic  materials ^  it  is  seen  that  the  displacement 
can  be  expressed  as 

y(x.t)  =  yttp — ^o^^^  , 

'  cr 

or  (1-6) 

«(*)  =  TT-p/p-- ■  «o  ' 

'  cr 

where  P  is  given  by  equation  (1.3)* 

O  L 

The  function  f(P,t)  represents  the  time  response  of 
the  viscoelastic  material.  If  for  convenience  we  assume 
the  load  is  applied  at  t  =  0  ,  then 

f(P,o)  =  1  ,  (1.7) 

and  v;e  have  from  equation  (1.6) 


y(x,o)  =  Y3p7p — y  (x)  ,  (1.8) 

'  cr 

which  is  identical  to  the  results  for  an  elastic  column  anal¬ 
ysis  [l],  using  the  same  assumptions.  It  is  obvious  that  as 
P  P  the  displacement  increases  very  rapidly,  and  thus 
P  is  an  upper  limit  on  the  loads  that  may  be  supported  for 
a  very  short  periot?  of  time. 


The  character  of  the  function  f(P,t)  depends  very 

P 


strongly  on  the  load  P  .  Let  P 


cr 


be  defined  as  the  long 


time  critical  load,  which  is  given  by 

P 


p  TT  EpI 


cr 


(1.9) 
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where  Ep  is  the  long  time  or  final  modulus  of  the  viscoe¬ 
lastic  material,  and  is  defined  as  the  modulus  governing  the 
stress  for  an  infinitely  long  time  application  of  a  constant 
strain.  A  viscoelastic  fluid  is  characterized  by  having 
=  0  while  for  a  viscoelastic  solid  0  <  E^  <  E(0)  . 

These  properties  are  more  fully  described  in  Chapter  2. 

For  „ 

P  >  Per  '  as  t-»  ,  (1.10) 

and  for 

p  l-P/Por 

p  <  p„/  ,  f(p,t)-. - — ^  as  t-.».  (1.11) 

1-P/Pcr 

P 

Thus  it  is  seen  that  for  loads  P  >  ,  the  column  is 

essentially  unstable  as  the  deflection  keeps  increasing  in- 

P 

definitely,  while  for  loads  P  <  P  the  column  is  stable 

in  the  sense  that  the  deflection  remains  finite.  Pig.  1.2b 

shows  a  qualitative  plot  of  the  column  deflection  vs.  time. 

P 

The  curve  for  P  <  P  is  asymptotic  to  the  dashed  line 

w  X 

=  - 1 - -  while  for  P  ^  <  P  <  P  „  the  dashed  line 

1-P/P^/  ' 

has  no  significance  whatsoever.  It  is  apparent  from  Fig.  1.2b 
that  loads  in  the  range  P^^‘  <  P  P^p  supoorted 

for  a  short  time,  but  for  long  time  loading  it  is  necessary 
that  P  <  P^^^  . 

Comparing  then  the  results  of  the  elastic  and  visco¬ 
elastic  analysis,  it  is  seen  that  the  critical  buckling  load 
of  a  perfect  column,  P  ,  serves  as  an  upper  limit  for 

V  JL 


-9- 


the  safe  load  thai:  can  be  supported  by  an  elastic  column, 
even  if  the  column  is  not  perfect.  However,  for  the  visco¬ 
elastic  column,  P  serves  as  an  upper  limit  only  for  the 

cr 

perfect  column  and  for  very  short  duration  loadings  of  the 
imperfect  column.  If  the  load  is  to  be  applied  for  a  long 
time  to  a  viscoelastic  column,  then  even  for  an  infinitesi¬ 
mally  small  imperfection  it  is  necessary  to  restrict  the  load 

ri 

to  be  less  than 

cr 

A  major  problem  in  describing  the  behavior  of  an  imper¬ 
fect  viscoelastic  column  occurs  because  in  the  analysis  there 
is  no  natural  criterion  with  which  to  associate  failure  other 
than  the  fact  that  for  some  load  conditions  the  displacements 
becane  infinite.  However,  the  displacements  approach  infin¬ 
ity  only  as  the  time  also  approaches  infinity,  and  so  this 
does  not  provide  a  unique  measure  of  the  column  response  for 
different  viscoelastic  materials.  Thus,  in  the  sense  of 
some  natural  criterion,  a  finite  critical  buckling  or  failure 
time  does  not  exist  for  any  value  of  the  load  below  P 

C  I* 

W,  [8]. 

The  second  kind  of  buckling  to  be  considered  is  of  the 
snap-through  type  as  exhibited  by  a  shallow  pin  ended  arch 
of  uniform  cross-section,  loaded  as  shown  in  Pig.  1.3a« 
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(a)  (b) 

Pig.  1.3 

In  discussing  the  response  of  such  an  arch  it  Is  convenient 
to  consider  a  loading  consisting  of  only  a  concentrated  load 
P  as  shown,  although  all  the  remarks  carry  through  for  dis¬ 
tributed  loads  as  v/ell.  The  deflection  of  zhe  midpoint  of 
the  arch  under  the  load  is  designated  by  6(t)  ,  thus 


where 

and 


^(t)  = 


=  , 


(1.12) 


defines  the  initial  unloaded  shape  of  the  arch. 


■^o  =  5^0  (?) 


(1-13) 


The  load-deflection  response  of  an  elastic  arch  [11] 
is  given  in  Pig.  1.3b.  Curve  A  corresponds  to  an  arch 
whose  initial  rise-to-span  ratio,  ,  is  such  that  its 

load  deflection  curve  has  a  horizontal  tangent  at  the  point 
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6/6^  =  1  ,  while  curve  B  represents  an  arch  whose  ini¬ 

tial  rise-to-span  ratio  is  less  than  that  of  the  arch  repre¬ 
sented  by  A  .  For  these  two  cases  the  load  deflection  re¬ 
lation  has  a  one-to-one  correspondence  and  there  is  no  in¬ 
stability  for  any  magnitude  of  the  load.  A  physical  expla¬ 
nation  for  this  behavior  is  that  because  the  arches  are  so 
shallow,  the  resistance  to  deformation  by  bending  is  much 
greater  than  that  due  to  the  axial  thrust,  and  so  the  loss 
of  the  resistance  by  the  axial  thrust  as  the  deflection  in¬ 
creases  is  more  than  offset  by  the  increase  due  to  bending. 

Curves  C  and  D  represent  arches  which  ar’e  perfectly 
symmetrical  with  respect  to  the  line  ^  ^  ,  and  whose 

rise-to-span  ratios  are  greater  than  those  of  A  .  In  these 
cases  there  is  a  range  of  loading  where  the  deflection  is 
a  multivalued  function  of  the  load,  and  so  there  is  more 
than  one  equilibrium  position  for  every  value  of  P  .  For 
chese  cases  then  there  is  a  well  defined  buckling  load  as 
shown  by  (P  )„  and  (P-^)-n  lob,  for  at  these 

loads  the  respective  arches  undergo  a  finite  deformation, 
as  shown  by  the  dashed  arrows,  to  a  new  stable  equilibrium 
position. 

Curve  C  represents  an  arch  whose  initial  rise-to-span 
ratio  is  such  that  when  the  arch  buckles,  the  deformations 
in  going  from  the  unbuckled  to  the  buckled  state  are  symmet¬ 
rical,  as  shown  in  Fig.  1.4a.  Curve  D  then  represents  an 
arch  whose  rise-to-span  ratio  is  greater  than  C  ,  and 
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which  buckles  unsymmetrically  as  shovm  in  Pig.  1.4b. 


Fig.  1.4 


The  term  snap-through  can  readily  be  seen  to  be  a  good 
description  of  this  type  of  behavior,  as  the  structure 
under  a  dead  load  would  physically  snap  through  to  the 
post-buckled  position. 

For  the  case  depicted  by  curve  C  the  equilibrium 
deformation  path  (curve  C)  is  seen  to  be  a  continuous 
curve  with  no  bifurcation  points.  However,  the  equilibrium 
deformation  path  for  case  D  consists  of  the  curve  D 
plus  the  dotted  portion  as  shown  in  Pig.  1.3b.  The  dotted 
portion  corresponds  to  symmetrical  deformations,  which  are 
equilitrlum  deformations  but  are  unstable.  Thus  at 
P  =  there  is  a  bifurcation  point,  just  as  there  was 

with  the  perfect  column. 

The  load  deformation  curves  for  imperfect  arches,  that 
is,  arches  whose  loading  or  initial  shape  is  not  symmetrical 
with  respect  to  the  line  x  =  ^  ,  are  shovm  in  Fig.  I.5 
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with  the  dotted  lines  the  corresponding  curves  for  the  sym¬ 
metrical  arches.  (Por  discussion  purposes  here  it  is  assum¬ 
ed  that  the  non-symmetry  is  caused  by  the  initial  shape,  and 
that  this  non-symmetry  could  be  expressed  by  the  second  mode 
of  a  Fourier  sine  expansion  of  the  initial  arch  shape.) 


Fig*  1*5 


The  equilibrium  deformation  paths  as  shown  are  now  all 
continuous  with  no  bifurcation  points,  but  there  still  exist 
well  defined  buckling  loads  ^^cr^D  ’ 

in  contra.st  to  the  imperfect  elastic  column  which  has  no 
well  defined  buckling  load.  Following  the  terminology  in 
[17],  unsymmetrical  buckling  at  bifurcation  points  will  be 
termed  transitional  buckling,  in  order  to  distinguish  it  from 
unsymmetrical  buckling  where  there  are  no  bifurcation  points. 
Buckling  where  no  bifurcation  point  is  present  will  be  termed 
non-transitional . 
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In  this  thesis  a  solution  is  given  of  the  viscoelastic 
arch  problem.  It  is  found  that  the  response  of  a  viscoe¬ 
lastic  arch  is  in  many  ways  s^.milar  to  an  elastic  arch.  If 
the  arch  is  very  shallow  then  there  is  no  instability; 
higher  arches  may  buckle  in  either  symmetrical  or  unsymmet- 
rical  modes j  depending  on  the  loading  and  the  initial  rise- 
to-span  ratio.  The  main  effect  of  the  viscoelastic  material 
is  to  permit  the  deformation  to  increase  with  time,  even  for 
a  constant  load  P  .  This  essentially  decreases  the  stiff¬ 
ness  of  the  arch  in  resisting  additional  load  increments, 
and  consequently  decreases  the  buckling  load.  Thus  it  is 
possible  for  a  constant  load  P  to  be  stable  at  one  time, 
but  at  some  later  time  unstable. 

In  contrast  to  the  response  of  a  viscoelastic  column, 
the  viscoelastic  arch  buckles  at  a  finite  time  if  it  buckles 
at  all.  For  a  sufficiently  small  load  and  a  viscoelastic 
solid  material  the  arch  may  always  be  stable,  or  it  may  be 
so  shallow  to  begin  with  that  instability  does  not  occur. 

It  is  this  latter  case  that  most  closely  resembles  the  col¬ 
umn  response,  since  the  arch  deflections  may  keep  on  in¬ 
creasing  indefinitely  v/ithout  buckling  action  if  the  visco¬ 
elastic  material  is  of  the  fluid  type.  Another  difference 
between  the  viscoelastic  arch  and  column  response  is  that 
for  an  imperfect  arch  whose  non-symmetry  is  of  infinitesimal 
magnitude,  the  buckling  load  and  buckling  time  is  different 
from  the  buckling  load  and  buckling  time  of  the  corresponding 
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perfect  arch  only  by  an  infinitesimal  amount.  It  was  seen 

Tp 

earlier  that  for  P  >  P  ,  the  response  of  an  imperfect 

cx 

viscoelastic  column  varied  widely  from  the  perfect  column 
no  matter  what  the  magnitude  of  the  imperfection. 

The  analysis  as  given  in  the  main  text  is  valid  for  an 
arbitrary  initial  arch  shape  (provided  the  arch  is  shallow), 
arbitrary  lateral  loading,  and  any  linear  viscoelastic  mater¬ 
ial.  Because  numerical  methods  are  required  in  the  solution 
of  the  governing  coupled  non-linear  integral  equations,  it 
is  possible  to  use  directly  measured  material  properties  in 
their  numerical  form. 

The  buckling  criterion  used  in  the  solution  is  based  on 
stability  with  respect  to  infinitesimal  displacements  about 
an  equilibriu  position.  At  any  instant  of  time  the  response 
to  an  instantaneous  displacement  is  governed  by  the  initial 
modulus  of  the  material,  and  so  the  stability  can  be  investi¬ 
gated  in  the  same  general  manner  as  the  stability  of  an 
elastic  arch,  provided  the  time  history  of  the  deformation 
up  to  that  instant  is  known.  The  theory  also  takes  into 
account  the  possibility  of  buckling  by  the  two  different 
mechanisms,  the  transitional  or  non-trans:' tional  modes  of 
failure . 

The  constitutive  relations  of  the  linear  viscoelastic 
material  are  expressed  in  integral  operator  form  (see  Chap¬ 
ter  2)  by  the  use  of  the  relaxation  modulus  E(t) 
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as  opposed  to  the  dif¬ 


ferential  operator  approach  ^P(a(t))  =  Q{e(t)) 

P  and  Q  are  linear  differential  operators  in 
the  arch  problem,  which  is  geometrically  non-linear,  tht- 
integral  operator  representation  of  the  material  offers 
several  advantages  over  the  more  common  differential  operator 
representation.  These  advantages  do  not  hold  for  the  visco¬ 
elastic  column  problem  since  the  column  problem  has  a  govern¬ 
ing  equation  which  is  linear.  The  most  widely  accepted  method 
of  solving  linear  viscoelastic  problems  is  by  the  use  of 
the  Laplace  transform  [13 and  in  the  transform  problem 
either  representation  of  the  material  properties  gives  ex¬ 
actly  the  same  result. 

Some  of  the  advantages  of  the  use  of  integral  operator 
relations  in  the  arch  problem  and  in  general  are: 

(a)  the  relaxation  modulus  can  be  measured  directly 
by  experiment. 

(b)  the  general  form  of  the  governing  equations  remains 
the  same  as  for  the  elastic  arch  problem.  This  allows  the 
same  type  of  stability  analysis  and  affords  a  good  physical 
understanding  of  the  viscoelastic  effects. 

(c)  the  numerical  integration  is  simplified.  This  is 
especially  true  if  the  material  is  quite  complex,  and  would 


,  where 
t  ^  .  For 
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require  the  use  of  high  order  differential  operators  if  this 
approach  were  used. 

Curves  showing  deformation  vs.  time  for  various  load 
levels  and  buckling  time  vs.  load  are  presented  for  two  dif¬ 
ferent  materials.  No  attempt  was  made  to  give  results  for  a 
wide  number  of  situations;  the  examples  were  chosen  so  as  to 
show  the  effect  of  some  asymmetry  in  the  structure  or  its 
loading,  and  to  illustrate  the  basic  difference  in  response 
between  viscoelastic  solid  and  fluid  materials. 

Arches  made  of  viscoelastic  solid  materials  may,  if  the 
load  is  small  enough  and  constant  as  t  ,  approach  a  final 
stable  equilibrium  position.  The  equations  governing  the 
equilibrium  position  in  this  case  reduce  to  the  elastic 
equations  with  a  modulus  of  elasticity  of  Ep  ;  these  are 
algebraic  equations  and  can  be  solved  to  find  the  final  e- 
quilibrium  position.  The  value  of  the  final  equilibrium 
position  can  then  be  used  to  evaluate  the  integrals  in  the 
equations  that  govern  the  buckling  action,  and  this  reduces 
these  equations  to  algebraic  equations.  It  is  then  possible 


to  determine  the 


j  j  u  ^  *1 

auuxuxuliax 


load  that  would  cause  buckling 


at  this  time.  Unfortunately  the  algebraic  equations  may  be 
of  quite  high  order  and  even  for  the  simplest  case  the  ad¬ 
ditional  load  cannot  be  expressed  in  explicit  form.  However, 
this  method  can  be  used  to  determine  the  Instantaneous  buck¬ 
ling  load  of  an  arch  which  has  reached  a  final  position  with¬ 
out  needing  to  know  the  entire  relaxation  function  E(t)  of 
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the  material  or  the  time  history  of  deformation.  All  that 
is  required  is  the  instantaneous  modulus  E(0)  and  the  final 
modulus  Ep  ,  'Which  could  actually  be  determined  from  the 
final  position  of  the  arch  if  this  viere  known. 

Results  for  this  long  time  problem  for  several  values 
of  the  final  modulus  and  various  loading  conditions  are  given 
to  compare  the  long  time  buckling  load  to  the  instantaneous 
or  elastic  buckling  load. 

Considerable  v;ork  has  been  done  on  the  so-called  creep 
buckling  of  columns.  Creep  buckling  is  usually  concerned  with 
materials  whose  constitutive  relations  are  non-linear,  as 
opposed  to  the  linear  constitutive  equations  used  in  linear 
viscoelasticity.  Most  existing  creep  laws  do  not  represent 
instantaneous  and  retarded  elasticity  as  v/ell  as  viscous  flow 
v/hich  are  all  properties  of  linear  viscoelasticity;  however, 
for  metals  it  is  generally  agreed  that  the  non-linear  creep 
lav;s  give  a  better  representation  of  the  material  behavior. 
Hoff  [l6]  gives  a  survey  of  the  theories  of  creep  buckling 
up  to  1958. 

Pian  [10]  has  given  an  analysis  for  the  creep  buckling 
of  a  uniform.ly  loaded  arch  using  a  2-element  non-linear  Kel¬ 
vin  model  to  represent  the  material.  By  linearizing  the  model 
his  problem  reduces  to  the  linear  viscoelastic  case.  The 
Kelvin  model  is  not  a  particularly  good  representation  of 
material  behavior,  but  because  of  its  simplicity  Pian  is  able 
to  get  a  numerical  solution  for  both  the  linear  and  non-linear 
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cases.  However,  his  solution  does  not  take  into  account  the 
possibility  of  transitional  buckling,  and  for  cases  where  tne 


transitional  mode 
ling  time,  Pian*s 


gives  the  critical  buckling  load  or  buck- 
results  predict  larger  values. 
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Ch-jpLer  II  LIr.«-  i-  Viscoelastic  Ma^ex‘Inic 

Viscoelast  ic  materials  are  char- ct ''r^  ;  ,7  th  l*'  pi-op- 

erty  to  deform  with  time  under  the  action  of  -  given  ioodlr.g. 
In  a  linear  viscoelastic  material,  in  particular,  the  con¬ 
stitutive  relation  is  such  that  the  ratio  of  stress  strain 
at  any  Inst'e  t  ‘s  ■  constant  IrKh.'pendent  of  the  raagidtude  ol" 
the  stress  and  strain. 

In  general  tv/o  independent  material  functions  ar'c  c-.- 
quired  to  define-  completely  a  homogeneous,  isotropic  visco¬ 
elastic  material  in  the  same  manner  that  two  independent 
materlaJ  const-'.nts  are  required  to  define  a  coi’respondii^g 
elastic  m-itvi'lai.  However,  hi  the  work  th- 1  i’oL  lows  only 
one  material  function  need  be  considered  because  of  the  one- 
dimensional  nature  of  the  problem.  For  more  general  visco¬ 
elastic  analysis  see  Bland  [12]  or  Lee  [I3]. 

For  this  problem  the  most  convenient  formulation  of  the 
stress-strain  relation  is  by  the  use  of  the  relaxation  mod¬ 
ulus  j£(t)  ,  which  Is  der'ined  by  the  relation 

E(t)  =  ,  (  .  I  ) 

where  a(t)  is  the  stress  resulting  from  an  applied  strain 

e(t)  =  t  (2.2) 

,  0,  t  <  O 

is  the  unit  step  function  H(t)  = 

1,  t  >0 


Since  the  material  is  considered  linear  this  relation  musr 

hold  for  any  value  of  e  .  E(t)  =  0  for  t  <  0  since 

o 

the  material  is  assumed  to  be  undisturbed  for  t  <  0  . 

The  general  form  of  E(t)  is  the  same  for  all  materials. 
There  is  an  initial  discontinuity  at  t  =  0  ,  and  the  func¬ 
tion  then  decreases  monotonically  to  some  final  value,  as 
shown  in  Pig.  2.1. 


Pig. 


C, 


1 


The  "initial  modulus"  of  the  material  is  defined  as 
where  =  E(0'^)  Similarly  the  "final  modulus"  is  de¬ 
fined  as  where  E^  =  11m  Eft^ 

p  P  ^  ^  — 

Linear  viscoelastic  materials  fall  into  two  general 
classifications,  fluids  and  solids.  Por  a  fluid  Ep  =  0  , 

whereas  for  a  solid  Ep  >  0  .  A  simple  model  representation 

1  + 

A  variable  written  as  t  refers  to  the  value  of  this 
variable  at  t  +  €  ,  and  f  =  t  -  e  ,  where  e  is  a 

positive  infinitesimal  quantity. 
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of  a  fluid  is  the  Maxwell  model.  Fig.  2.2,  in  which  a  lin¬ 
ear  dashpot  is  connected  in  series  to  a  linear  spring. 


K 


c(t) 


A/\A 


1  +  6{t) 


a(l) 


Pig.  2.2  Maxv^fell  Fluid  Model 


fluid 


Let  t  represent  real  time,  then  for  the  Maxwell 

E 

~ir  ^ 

E(t)  =  E^e  ^ 


(2.3) 


Define  a  dimensionless  time  parameter  t  as 

E 


(2.4) 


and  redefine  the  relaxation  modjlus  as 


E(t)  =  e""'  .  (2.^>) 

E 

“  is  termed  the  relaxation  time  of  the  Maxwell  material. 

The  simplest  model  representation  of  a  solid  material 
that  exhibits  an  initial  modulus  is  the  three  parameter  solid 
of  Fig.  2.3* 
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For  this  model 


ihllo) 


E(t)  =  E^e 


?r— ^  t  SiE„  / 

”  -e^Ve-^-^ 

1  O  \ 


(h  ^ 


.  (2.6) 


E,  +  E 

The  characteristic  relaxation  time  is  - - — -  and  so 

let  E  +  E 

T  =  _i - S  t  . 


(2.7) 


The  final  modulus  is  clearly 


E  E 

Ep  =  lim  E(t)  =  2  +  E 

t  CO  1  o 


(2.8) 


In  its  most  convenient  form  the  relaxation  modulus  is 
then  v;ritten  as 

E(t)  =  E^  e“  ■"  +  Ep(l  -  e"  '^)  .  (2.9) 

The  essential  difference  in  the  model  for  a  fluid  or 
solid  is  that  for  a  fluid  there  is  at  least  one  dashpot  con¬ 
nected  in  series,  whereas  the  dashpots  in  a  solid  material 
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I 


are  always  connected  in  parallel  with  a  spring. 

In  the  work  that  follows  the  time  parameter  t  is  al¬ 
ways  used.  For  real  materials  with  measured  relaxation  prop¬ 
erties,  it  is  necessary  to  choose  a  characteristic  relaxation 
time  in  order  to  define  a  dimensionless  time  parameter,  but 
this  can  always  be  done.  Lee  and  Rogers  [l4]  discuss  the 
use  of  measured  material  functions  in  stress  analysis  prob¬ 
lems  . 

Since  the  stress-strain  relations  are  linear  the  prin¬ 
ciple  of  superposition  is  valid  and  the  stress  for  a  smooth¬ 
ly  varying  strain  can  be  represented  by  the  integral 

T 

a(T)  =  U  ^  e(|)d^  .  (2.10) 

»  00 

Such  an  integral  is  known  as  a  Duhamel  integral,  and  the 
kernel  E(t)  is  related  to  the  Green's  function  of  the  dif¬ 
ferential  operators  that  can  alternately  be  used  to  express 
the  stress-strain  relations. 
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Chapter  III  General  Analysis 


Consider  a  shallow  pin-ended  arch  of  uniform  cross- 
section  loaded  as  shown  in  Fig.  3*1*  The  lateral  load 
q(x,T)  always  remains  iirected  parallel  to  the  y  axis. 
Before  the  application  of  the  lateral  load  the  arch  is  con- 
lidered  to  be  unstrained  and  the  location  of  its  centroidal 
axis  is  given  by  YqCx)  .  Under  the  action  of  the  lateral 
load  the  centroidal  axis  v/ill  be  displaced  to  a  new  position, 
and  let  this  be  denoted  by  y(x,T)  .  Without  loss  of  gen¬ 
erality  we  can  assume  q(x,T)  =  0  for  x  <  0  ,  and  so 


y(x,0-)  = 


(3.1) 


becomes  the  initial  condition  for  y(x,x)  .  P(t)  is  the 
axial  force  induced  in  the  arch  by  the  end  reactions  and 
is  directed  as  shown. 

Assuming  that  [y^j  and  |yj  are  much  smaller  than 
L  ,  that  the  curvature  remains  small  at  all  times  so  that 

is  negligible  in  comparison  to  1,  and  that  the 

thickness  of  the  arch  is  much  smaller  than  its  radius  of 


curvature;  then  the  usual  beam  theory  gives 


€(x,C.t)  = 


^  y(x,T)  d  y^(x) 


^  ,  (3.2) 


where  €(x,Cjx)  is  the  strain  of  a  longitudinal  fibre  of 
the  arch  located  the  distance  C  from  the  centroidal  axis 


as  shown  in  Pig.  5.2,  and  ^^(t)  is  the  strain  of  the  oen- 

troidal  axis  v;hich  is  assumed  constant  along  the  length  of 
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the  arch.  This  assumption  relies  on  the  fact  that  we  are 
considering  a  shallov/  arch  and  that  we  also  assume  the  axial 
force  in  the  arch  to  be  independent  of  x  .  The  moment 
M(x,t)  as  shown  in  Pig.  3*2  is  given  by 

-  iM(x,t)  =  JJ  C  a(x,C>T)dA  ,  (3-3) 

k 


where  the  negative  sign  is  introduced  to  set  up  tne  sign 
convention  that  positive  moment  causes  compressive  stresses 
on  the  positive  C  side  of  the  arch. 

Now  introduce  the  constitutive  relation  based  on  equation 

(2.10)  ^ 

=  J e(t  -  t)  §1  6(x,C<5)'i4  .  (3-'t) 

0“ 

and  substitute  this  into  equation  (3*3) ^  which  gives 

T 

-M(x,t)  =  JJ  C,  J E(t  -  4)  ^  e(x,C.4)dC  dA  .  (3.5) 

A  o~ 


Interchanging  the  order  of  integration  and  using  equation 
(3-2)  results  in  the  moment-curvature  relation 

T 

+  M(x,t)  =  ij  E(t  -  i) 
o~ 


where  I 


The  moment-load  relationship,  using  simple  arch  theory, 
is 
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(3-7) 


=  .  _^('p(^)y(x,x))  -  g(x,T) 

ax'^  ax'-  ^  ' 

for  the  loads  as  directed  in  Pig.  3*1> 

Since  the  arch  supports  are  assumed  rigid  the  projected 
length  L  of  the  arch  must  remain  constant,  thus 


0  =  - 


+ 


where  the  integral  in  the  above  expression  is  zo  a  first  ap¬ 
proximation  the  change  in  projected  length  caused  by  the  de¬ 
flection. 

The  axial  force  p(t)  is  given  by 

P(t)  =  -  JJ  cy(x,C>T)dA 

A 

and  by  using  equation  (3-^)  this  becomes 


I 

P(t)  =  -  -  ?)  ^  €(x,C,5)ci5  dA 


A  o 


If  the  order  of  integration  is  interchanged  and  equation  (3* 2) 
is  substituted  into  the  resulting  equation,  it  then  yields 


P(.)  =  -/e(x  [s(t)  - 


dA  d^ 


T 

=  -  a/e(x  -4)  ^  ,  (3.9) 

o" 
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1 


since 


A=/'/' 


JJ 
A 


dA  .  Now  by  using  equation  (3*8)  we  can  ex¬ 


press  P(t)  in  terms  of  the  position  of  the  arch;  thus 

•  (3-10) 


If  equation  (3-6)  is  differentiated  twice  with  respect 
to  X  and  substituted  into  equation  (3*7)^  we  have 


/e( 


T  - 


d^  =  -  P(t)  -  -  q(x,T)  . 

^x'" 


By  using  equation  (3* 10),  P(t)  in  the  above  equation  can 
be  eliminated  which  results  in 

_  A  ^^y(x,T) 


2L 


ax 


+  iJe{-v  -  i)  =  -  q(x,T)  ,  (3-11) 

the  governing  incegro-differential  equation  for  the  arch. 

The  boundary  conditions  are 


y(o,T) ■=  y(L,T)  =  0 


p  y(ojT)  _  ^^y(L,T)  _ 


(3.12) 


3x^ 


ax‘ 


=  0 


and  the  initial  condition  is 


y(x,o-)  =  y^Cx)  . 


(3.13) 


-3C- 


If  and  y(x,T)  are  expressed  in  a  Fourier  sine 
series,  this  will  satisfy  the  boundary  equations  (3.12),  and 
allow  the  space  derivatives  in  equation  (3.11)  to  be  removed. 


Thus  let 


y  (x)  =  2  a  L  sin 


(3.1't) 


y{x,T)  =  S  b^(T)L  sin  2=^  , 


(3.15) 


v/here  a  is  a  dimensionless  constant  and  b  (t)  is  a  di- 
mensionless  function  of  time. 

It  is  convenient  also  to  express  the  lateral  load 
q(x,T)  in  a  Fourier  sine  series;  thus 


q(x,T)  =  q^  2  ^ 

m=l 


{3.16) 


where  R  (t)  is  a  dimensionless  function  of  time  and  q^ 
is  a  constant  having  the  same  dimensions  as  q(x,T) 

Substituting  equations  (3.1^) ^  (3.15)  and  (3.l6)  into 
equation  (3.11 )  gives 

S  A„(x)  slnSS\fE(t  -  S  lcV(4)}d5 


/e{. 


d  fw  ~  4  imrxl  ,> 

^  173  "i”’  — r 

m=i 


m=l 


(3.17) 


-31- 


where 


/H 


-f( 


A2 


2  mTrb_(T)  cos  21|!2Sj  (jx 
\n=l  “■  -  / 


_2t  00 


F  L  “  ,2,  2,  X 


because  of  the  orthogonality  property  of  the  cosine  functions. 

Interchanging  the  order  of  integration  and  summation  in 
the  second  term  of  equation  (3-17)>  and  rearranging  results  in 

4, 


CO 


2  sin 
in=l 


mine 


+ 


J  1^(t  -  4)  ^  \{i)  de  +  qo\(T)  1  =  0 


Because  the  sine  functions  are  orthogonal  over  the  interval 


o  to  L  each  coefficient  of  sin 
vanish,  and  so 


mirx 


in  the  series  must 


.2  Y  A  K  r 

‘  JlT  J 


t\  A. 
d? 


CO  2  2 

-  V-/ 


.k=l 


m 


4  F^I 


L- 


/e{t  -  e)  ^  b|^(e)  d5  =  -  qgR^(t)  (3.18) 


m  =  1 ,  2,  3  • •  •  • 

It  i.s  now  convenient  to  introduce  new  notation  to  sim¬ 


plify  the  above  equations.  Let 


* 


A 


m 


-k  /K 


d^-  i  2^1 


,  and 


^o 


4 

2Tr  E  I 
o 


'1 

A 


Substituting  these  into  equations  (3.I8)  results  in 


-  e)  _d  “  ,.2^  2 


,,  z  k"B  "(e)de 

k=l  ^ 


R  (t) 
m  ^  ^ 

2 

m 


m  =  1,  2,  3  .... 


(3.19) 


with  the  initial  conditions 

B^(o  )  —  Aj^  n»  =  i,2^3***»  (b*20) 

Here  E(t)  ,  R  (t)  and  A_  are  known  in  the  problem 

m  m 

and  it  remains  to  find  determine  the  stability 

of  the  arch. 

The  loading  funcx^ion  R  (ti  will  normally  be  dlscon- 

m 

sinuous  at  t  =  0  (since  R  (t)  =  C  for  t  <  0  ),  and 

since  we  are  assuming  quasi-static  response  will  also 

be  discontinuous  at  t  =  0  as  long  as  E  .  It  is  con- 

o 

venient  to  remove  this  discontinuity  at  t  =  0  by  performing 
the  integration  in  two  steps,  from  0"  to  0+  and  0+  to 
T  .  Thus  equations  (3-19)  become 


m 


0+  ° 


R  (t) 

m 

2 

m 


m  =  1,  2,  3 


Integrating  the  above  integrals  by  parts,  and  by  using 
the  initial  conditions  as  given  by  equations  (3*2;0),  yields 


m  -  1 ,  2,  3  • • • •  (3*21) 


where  the  prime 
argument,  i.e. 

At  T  =  0+ 


denotes  differentiation  with  respect  to  the 

e'(t  -  o  =  ‘'div  -r-||  • 

,  equations  (3-21)  reduce  to 
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R  (0+) 

2 

m 


m  A 


m 


) 


m  =  1 ,  2,  j  .  .  . 


(3.22) 


Which  correspond  exactly  to  the  equations  of  the  purely  e- 
lastic  arch  of  modulus  [11]. 

Equations  (3-21)  are  an  infinite  system  of  coupled  non¬ 
linear  integral  equations,  and  as  such  give  very  little  en¬ 
couragement  for  finding  a  closed-form  solution.  However,  it 
is  possible  to  integrate  these  equations  numerically,  since 
it  will  be  shown  that  as  long  as  the  arch  is  stable  any  mode 

that  is  not  excited,  i.e.  A  =  R  (t)  =  j  ,  will  have  only 

m  m 

the  trivial  solution  B^(t)  =  0  .  Thus,  as  long  as  only  a 

finite  number  of  modes  are  excited,  our  system  can  be  reduced 
to  considering  only  this  finite  number  plus  one  more. 

To  integrate  equations  (3.21)  numerically,  replace  the 
time  integrals  with  a  step-by-step  sun....ation.  Therefore, 
let  ~  ^  ('^Q  =  O’*")  where  n  =  t.he  number  of  time  steps 

required  for  the  time  to  progress  from  0"^  to  t  . 

Then 


n— 1  ^i  +  1 

n— 2  p  i.  +  1 

2  / 

=  ^  / 

ir:0 

i=c  ^ 

•^i 

'^n-l 

(3.23) 


By  using  a  simple  trapezoidal  approximation,  the  convolution 
type  integrals  v/e  are  concerned  with  can  be  expressed  in 
the  following  form: 


f 


d 

‘1 


i+i  E‘(t  -%) 

f(0ci4 


f(T^)  +  r{T^^^) 

L  ^o 

2 

.  «. 

(5.24) 


Applying  equations  (3.23)  and  (3.24)  to  the  first  inte¬ 
gral  in  equations  (3.21),  one  gets 


E_  \k=l 


E(t  - 
'  n 


-  2A, 


+  ^“2  E(t^-t^)  -  E(t^- 


i=0 


E 


(k=l 


(3.25) 


Where 


E  ( T  -  T  -  )  -  E 
Vl)  =  — ^ 


2E. 


(3.26) 


ana 


-*6- 


I 


11-2  E(t„-t,)  - 


-i-  s 
i=0 


n  1 


(vi  o  /  Cl 


(3.27) 


Similarly  the  second  integral  in  equations  (3 ‘21)  gives 

/"  -  A^)de  =  +  sjT^)  , 


13-28) 


where 


n-2  E(t  -  T. )  -  2(t  -  V 

+  Z  - - - - - - - (t.)  +  B^(t.  .,  ) -2a'\  . 

i=l  2E  \^m'  i'  m'  i+1^  ry 

(3.29) 

Thus,  by  substituting  equations  (^.25)  and  (3.29)  into  equa¬ 
tions  ('j.2l),  one  obtains 


■75 —  +  m'^A  -  ms  ('T  )  m  =  1,  2,  3  ...  (3 '30) 

2  m  m'  n  / 


For  n  =  0  (t  =CH-),  equations  (3 *30)  reduce  to  equations 

(3.22). 

Assuming  that  the  arch  has  been  stable  up  to  and  in¬ 
cluding  T  =  >  then  at  the  next  time  step,  ~  '^n 

S  (t  )  and  S(t^)  are  known,  since  they  contain  only  the 
m  n  n 

previously  determined  values  of  B  (t  )  .  Thus  B  (t^) 

m  n'  m  n 

can  be  determined  from  the  system  of  equations  (3.30),  and 
a  check  can  be  made  of  the  stability  of  the  arch  at  this 
time , 

For  demonstrative  purposes,  two  different  loading  cases 

will  be  considered.  The  first  is  a  sinusoidal  arch  with  a 

sinusoidal  load,  i.e.,  =  R  (t)  -  0  for  m  >  1  ,  so 

that  only  the  first  mode  is  excited;  the  second  case  will 

have  an  additional  mode  excited,  i.e.,  one  A  or  R_(t)  /  0 

m  m 

for  m  >  1  ,  in  addition  to  A^  and  Rj^(t)  /  0  .  The  sec¬ 

ond  case  can  then  be  generalized  to  include  any  number  of 
excited  modes . 
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Chapter  IV  Sinusoidal  Arch  under  Sinusoidal  Loading 

Consjder  the  simple  case  of  a  lew  sinusoidal  arch  sub- 
jecced  to  a  sinusoidal  load  distribution,  then 


y^(x)  =  a^L  sin  ^  2/^  sin 


TTX 

L  ^ 


(4.1) 


?Tr*E.I 


q(x,t)  =  qQR^(T)  sin  ^  sin  ~  , 


Equations  (5-30),  the  general  equations  of  equilibrium, 
then  reduce  to 

-  +s(t„)]  =  -  R^(t„)  -  S,(T„) 

B2(Tr,)[(l  +  +  1  kV(^n)) 

-  -  S(T„)]  =  -  4S2(.„)  f 


m  n 


(l  +  E*(T^,T^.,))(n.2  +  1  kV(^n>) 

V  1  P  ' 

-  +  S(t^)  =  -  m  where  m  =  2,  3,  4  ... 


and  for  n  =  C  these  then  become 
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B^(0+)  1  +  Z  -  A,“  j=  -  R^(0+)  + 

B.(0+)  4  +  Z  A,  ^(0+)  -  A,^  =0 
^  L  1  ^  i  J 

• 

B^(0+)  jin^  +  2  A^2{0+)  -  =  0  . 


Equations  (4.3)  correspond  exactly  to  the  elastic  arch  buck¬ 
line  problem  and  are  discussed  at  length  in  [ll].  One  so- 
lutio'^  of  the  system  v/ould  have 

BjCO+j  =  82(0+)  =  Bi,{0+)  =  ...  =0 

B^(0+)ri  +  B,''’(0+)  -  =  -  R^(0+)  +  Aj 


The  relation  betv/een  B-  ’0+)  and  R,  (0+)  as  given  by 

X 

equations  (4.4)  is  plotted  in  Fig.  4.1. 

Depending  on  the  value  of  A^  there  are  several  pos¬ 
sibilities.  If  A^  <  1  ,  tne  curve  has  a  monotoric  slope  and 
so  there  exists  a  unique  one-to-one  relation  between  B^(0+) 
and  R^(0+)  .  Thus  for  any  give  loading  there  is  only  one 
equilibrium  position  and  the  system  is  stable.  If  A^  >  1  , 
then  for  R^(0+)  between  the  points  b  and  c  there  exist 
three  possible  equilibrium  positions.  Assuming  that 
B^(0")  =  A^  the  stable  equilibrium  position  of  the  arch  is  gi 
en  oy  the  portion  ab  of  the  curve,  provided  R, (0+)  <  R  (0+) 

X  S 


If  R,  (0+)  >  R  (0+)  the  only  ecjuilibriurr  position  is  along 

X  S 

de  .  Thus  it  is  apparent  that,  if  the  load  is  increased  be¬ 
yond  R,(0+)  the  arch  snaps  through  from  b  to  d  and  fol- 
lows  the  de  portion  of  the  curve. 

The  portion  be  of  the  curve  represents  unstable  equi¬ 
librium  positions,  while  along  cd  the  positions  are  again 
stable;  however  since  the  arch  starts  at  a  the  portion  ab 
represents  the  stable  unbuckled  equilibrium  positions,  and 
hence  the  point  b  is  the  critical  point  at  which  buckling 
occurs . 


R  (O"^)  ,  which  is  termed  the  non-transitional  buckling 
s 

load,  can  be  determined  by  setting 


dRi(0+) 

dB^(0+)  "  ^ 


d^R, (0+) 

for  — o -  <  0 

dB.  (0+) 

X 


(^.5) 


^/hich  results  in 


Rg(0+)  =  A-  + 


A^  >  1 


(^.6) 


Now  a  second  so:..ution  of  one  system  of  equations  (^.5) 
could  exist.  If  one  3jj,(0+)  in  addition  to  B^(0+)  is  dif¬ 
ferent  from  zero,  then  the  only  two  equations  of  the  set  (4.3) 
not  identically  satisfied  are 


B3^(0+)  jl  +  B^^(0+)  +  -  A^^j  =  -  R^(0+)  +  A^ 

B^(0+)jm^  +  B^^(0+)  +  -  A^^j  = 


0 


This  solution  leads  to  the  transitional  buckling  load  if 

B  (0+)  /  0  is  found  to  be  a  valid  solution.  If  it  is  not 

a  valid  solution,  then  the  arch  will  buckle  in  a  symmetrical 

mode  and  R„(0+)  will  be  the  buckling  load, 
s 

Since  it  is  assumed  that  B  (0+)  /  0  the  second  of 

equations  (4.7)  gives 

m^  +  B^^(0+)  +  m^B^^(0+)  -  =  0  •  (4.8) 

At  this  point  it  is  seen  that  only  one  B^CO't)  in  ad¬ 
dition  to  B^(0+)  can  be  non-zero.  If  it  is  also  assumed 
that  B^(0+)  0  ,  n  /  l,m  then  it  is  necessary  that 

m^  +  Bt^(0+)  +  m^B  ^(0+)  +  n^B^^(0+)  -  A,^  =  0  , 

i  m  n  '  1 

and 

V?  +  Bt^(0+)  +  mS  ^(0+)  +  n^B  ^(0+)  -  A^  ^  =  0  . 

1  m  n  '  1 


This  clearly  is  a  c-^ntradiction  since  n  /  m  and  so  the 
assumption  that  B^(0+)  /  0  is  incorrect. 

Substituting  equation  (4.8)  into  the  first  equation  of 
(4.7)  results  in 


\\j  -  ] 


Rl(0+)  +  A^ 


m 


-  1 


and  then  by  using  this  result  in  equation  (4.8),  one  obtains 

2 


m^B  ‘^(0+)  =  -  +  k‘ 

m  1 


'R,  (0+)  -  A, 


r\ 

m.  -  1 


(4.10) 
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Equatloi.  (4.9)  is  plotted  in  Pig.  4.2a  as  the  ful]  line,  where 
the  dotted  line  is  the  same  curve  as  that  of  Fig.  4.1.  Equa¬ 
tion  (4.10)  is  plotted  in  Pig.  4.2b. 

Prom  Fig.  4.2b  or  equation  (4.10)  it  is  seen  that 
B^(0'^)  has  a  realj  non-zero  solution  only  in  a  definite  range 
of  -^^(0+)  >  and  for  this  solution  to  exist  6^(0+)  must  lie 
on  the  stra'’.£,ht  lino  segment  b^c‘  of  Pig.  4.2a.  For 
B  (C"^)  =  0  the  solution  for  B,  (0+)  lies  on  the  dotted  curve 
of  Pig.  4.2a.  Thus  for  R^(0+)  <  R^(0+)  ,  B^(0'^)  is  given 
by  the  curve  ab’  ^  and  B  (0+)  =  0  .  For  R^ (0+)  =  R  (0+) 
the  arch  will  snap  througn  to  the  point  d*  ,  v;ith  B.|  (0+) 

following  the  path  b*cM*  and  B  (0+)  following  either  of 

m 

the  closed  paths  fghkf  of  Fig.  4.2b.  When  B^(O'i)  is  on 

the  path  b^c’  ,  B  (0+)  is  non-zero  and  the  arch  deflection 

m 

is  unsymmetrical . 

Thus  we  conclude  that  R^(0+)  =  R^^(0+)  the  critical 
transitional  buckling  load.  The  lov/est  possible  value  of 
R^(0+)  v/ill  occur  v;hen  m  =  2  is  substituted  into  equation 
(4.9)^  and  so  the  second  mode  v/ill  be  the  other  non-zero  mode 
during  the  unsymmetrical  snapping-through  action. 

In  order  for  the  arch  to  buckle  unsymmetrically  or  as  we 
have  defined  it,  transitional  bucxling,  the  point  b’  musi. 
lie  on  the  portion  ab  of  the  dotted  curve  of  Fig.  ^.2a,  oMi -r- 
wise  b  will  be  the  critical  point  and  the  arch  will  buckle 
symmetrically.  It  is  possible  for  b’  to  lie  on  that  pore fen 
of  the  dotted  curve  th'^t  goes  from  b  to  the  R^  (0+)  axis. 
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If  this  is  the  case  then  R^(0+)  is  the  critical  load  but 
the  deformation  in  going  from  b  to  d  v;ill  be  symmetrical 
from  b  to  b’  ,  unsymmetrical  from  b‘  to  c^  and  s^/m- 
metrical  again  along  c’cd  .  Since  the  critical  load  is  not 
changed  from  the  symmetrical,  non-transitional  buckling  load, 
and  since  the  buckling  is  initiated  in  the  manner  of  symmetri¬ 
cal  buckling,  it  v/ill  be  considered  as  non-transitional  in 
the  remainder  of  this  dissertation. 

Prom  equations  (4.^0  (^-*9)  it  is  seen  that  if 

Aj  ^  y/  5*  5  ,•  b*  will  fall  on  the  portion  ab  and  the 
critical  load  will  be  (0+)  =  R  (0+)  .  R,  (0+)  is  deter- 

mined  by  setting  B^(0+)  =0  in  equation  (4.10‘1  v/hich  gives, 

for  m  =  2  _ 

R^(0+)  =  ^  (^-11) 


Thus  uhe  critical  value  of  rhe  loading  is: 


R^(0+) 


=  R  (0+)  =  A,  + 
cr  s'  1 


i.'-  2 
27 


(Af  -  D' 


R3^(0+) 


cr 


for  1  <  A^ 


f—T - 

=  R^^(0+)  =  A^  +  3  /  A^^  -  4 

for  A,  >y~5-5  • 


Nov;  assuming  that  R^(0+)  < 


Ri{o+)] 


cr 


(4.12) 


,  i.e.  the  arch 


is  stable  at  t  =  0+  ,  is  given  by  the  second  of 

equations  (4.4)  and  B,(0+)  =0,m=2,  3,  4... 

fii 
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Thus 


% 


\ 

f 


the  important  conclusion  is  reached  that,  if  the  arch  is 

stable,  the  sinusoidal  configuration  is  maintained. 

Proceeding  on  to  the  next  time  step  t  =  ,  it  is  nov/ 

possible  to  calculate  3(t, )  and  S  (t. )  since  the  values 

1  mi 

of  B  (0+)  are  known.  Hence  from  equations  (3«27)  and 
m 

(3.29)  it  is  seen  that  S(t^)  and  non-zero, 

but  that  =0  ,  m  =  2,  3^  ^  ^  since  B^(0+)  /  0  , 

B  (0+)  =0  for  m  =  2,  5^  ^  •••  • 

Inserting  these  results  in  equations  (4.2)  gives,  for 


=  -  Ri(tP  +  -  Sj^(tP 


('*•13) 


(3  +  E*(Tj^,0+)J(m2  -  Aj^2  +  S(Tj^) 


=  0 


m  =  2,  3}  ^  • • •  • 


Equations  (4.15)  can  only  be  used  to  determine  the  sta- 
b''e  equilibrium  values  of  B^(t^)  under  the  assumption  that 
the  arch  does  not  buckle  unsymmetrically  in  a  transi^;ional 
mode,  and  provided  tnat  an  unbuckled  equilibrium  position 
exists.  Using  the  same  argument  as  v;as  used  for  t  =  0+  , 
this  requires  that 


R^(t^)  < 
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If  does  not  satisfy  this  Inequality,  then  the 

arch  has  become  unstable  sometime  between  i:  =  0+  and 
T  =  ,  and  has  snapped  through  to  a  buckled  equilibrium 

position.  Hov^ever,  if  equation  (4.l4)  is  satisfied,  the  un¬ 
buckled  equilibrium  ^alue  of  is  the  largest  root  of 

equations  (4.13)  when  ~  ^  ’  m  =  2,  Z>  4  ...  . 

Instantaneous  changes  in  deformation  are  governed  by 
the  initial  modulus  ;  equations  (4.13)  are  not  applicable 
to  this  type  of  response  since  they  are  valid  only  for  a 
small  step  in  a  smoothly  varying  continuous  change,  and  there¬ 
fore  they  cannot  be  used  to  predict  the  stability  of  the  arch 
v/ith  respect  to  the  transitional  buckling  mode.  They  are  va3- 
id  in  the  non-transitional  case  since,  to  the  accuracy  of  the 
numerical  v/ork,  they  give  the  equilibrium  position  and  if  no 
nearby  equilibrium  position  exists,  then  the  arch  must  have 
snapped  through  to  its  buckled  position. 

To  determine  if  the  arcn  is  stable  at  t  =  ,  it  is 

then  necessary  to  look  at  the  instantaneous  response  of  the 
arch  at  this  time.  This  car  be  done  by  applying  an  additional 
load,  say  ,  at  t  =  and  then  investigating  the  sta¬ 
bility  at  T  =  T,  by  considering  the  B  vs. 

i  mi 

{R^(t^)  +  J  relations.  In  order  that  this  additional 
load  not  affect  conditions  at  a  later  time  it  must  be  applied 
for  an  infinitesimal  time  only,  and  in  this  respect  it  can  be 
considered  as  an  imaginary  or  pseudo  load  used  exclusively 
for  the  determination  ol  stability. 
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Thus  under  the  assumptions  that  the  arch  does  not  buckle 


unsymmelricaily  and  that  an  unbuckled  equllibrjum  position 
exists,  found  from  equations  (4.13).  It  is  then 

possible  to  calculate  S(t^  )  and  )  to  be  used  in 

investigating  the  stability  at  t  =  .  By  definition, 

see  equation  (3.28) 


hov/ever 

E*(ti  +  ,Ti)  =  - ^  ^  =  0 


2E. 


2E 


ana 


o  • 


for  the  integrands  considered  here.  Thus 

Tn 


0+  ° 

and  by  using  equation  (3*28)  again  this  can  be  written  as 

\ 


A  similar  argument  follows  for  S(T^'t)  and  so 


S(t^  +  )  =  S(t;j^)  +  E*(t^,0+) 


>(^.15) 


S  (t  +)  =  0 
m '  1 


m  =  2,  3> 


•  • 


Substituting  equations  (4.15)  directly  into  equations 
(4.2),  and  noting  the  load  at  t  =  t^+  is 
gives 

Ei(ti  +  )  1  +  S  A^^{tj^+)  -  A,®  +  S(t^  +  )J  ^ 

L.  /C~l 


■■  ■  ^i^'^i^  ^i) ^1  " 

I 

B„(Ti+)p  =  0  j 


)(4.16) 


m  =  2,  3^  4  . . . 


Equations  (3-30)  would  give  the  same  result  if  it  v/ere 

considered  that  =  t^+  ,  "^n  1  ^  ^  ’  They  have  not  been 

used  since  it  is  lesired  that  d®  a  finite  time 

step,  and  it  is  convenient  to  write  t  -  x  as  an  inf  ini- 

n  i* 

tesimal  step. 

The  investigation  of  equations  (4.l6)  follows  the  gen¬ 
eral  procedure  used  in  looking  at  equations  (4.3).  Again 
two  possible  solutions  exist  and  the  load  deflection  relations 
are  plotted  in  Pigs.  4.3a  and  4.3b. 
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Pig.  4.3 


The  critical  loads  are  given  by 


Ri(t^)+Ar3  =R,(-r^+) 

^  ^  O  X* 


-  -  S^(tj^  +  )  +  ^27^1  ~  ^ 


for  1  <  -  S(t3^+)  <5-5 


[r,  (t,  )  +  AR J  =  R  (t,  +  ) 

L  -  J-Jcr 


=  A  -  S^(t+)  +  -  S(t^+)  - 


for  a/  -  S(t  +)  >  5.3 


If  ^  then  the  assumption  that 

the  arch  does  not  buckle  unsymmetrlcally  is  incorrect,  and  the 
arch  has  become  unstable  sometime  between  t  =  0+  and  t  = 

Thus  if  R^(t^)  <  R^(t^)  +  ^2.  J  ^  assumption  made 

earlier  regarding  the  existence  of  an  unbuckled  stable  equi¬ 
librium  position  is  valid,  and  the  arch  is  stable  under  the 
load  R^(t)  up  to  and  including  x  =  ,  for  R^(x)  monctoni- 

cally  increasing.  If  R^(x)  is  a  decreasing  function  of  time, 
there  is  a  possibility  of  the  arch  becoming  unstable  between 
T  =  0+  and  T  =  x^  even  if  the  analysis  shov;s  it  to  be  sta¬ 
ble  at  the  two  end  times  0+  and  x^  .  Thus  special  care 
must  be  taken  in  the  choice  of  time  steps  to  insure  that  this 
does  not  occur.  Of  course  if  the  arch  were  stable  at  the  time 
x^+  for  a  loading  ^  then  this  would 

insure  stability  during  tne  interval  0+  to  x^  . 

Once  stability  is  insured  at  x  =  x^  it  is  possible  to 
compute  S(xo)  ^  S  (xp)  and  proceed  on  to  the  second  time 
step,  where  once  again  the  same  analysis  can  be  made  to  deter¬ 
mine  the  conditions  at  that  time.  This  process  can  be  repeated 
until  the  arch  becomes  unstable. 

Pig.  4.4  shows  the  load  deflection  relation  for  a  partic¬ 
ular  arch  at  severax  instants  of  time,  and  illustrates  very 
clearly  the  effect  of  the  creep  on  the  response  and  buckling 
load  of  the  arch  as  the  time  increases. 


Pig,  k.k  Maxwe31  Fluid  Material 

Sinusoidal  Arch  A^  =  4 


Sinusoidal  Load  R^(t)  -  8 
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Chapter  V  Non-slnusoldal  Arch 

(V.l)  Tvio  modes  excited 

In  the  preceding  chapter  only  the  first  mode 
excited^  and  it  v/as  shown  that  if  the  arch  was  stable,  all 
other  modes  remained  zero.  Now  consider  an  arch  which  has 
one  additional  mode  excited,  say,  the  r  mode.  This  re¬ 
quires  that  or  R^(t)  or  both  be  non-zero. 

Equations  (3 ‘30)  then  become 


2  2  ^ 

-  K  -  rA  “  + 
1  r 


■  k-1 


.  2  2.2,-, 
-  Al  -  r 


- 


ana 


>(5.1) 


-  A,^  -  r^A  ^  +  S(t  )[  =  -  m^S  (t  ) 

1  r  “J  mn 

where  m  l,r  . 


For  n  =  0  these  reduce  to 
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Bi(0+) 


Br(0+) 


“>  ?  P  /  .  \  2  22" 

1+2  k^B-  "^(0+)  ~  K  -  ^  = 

k=l  ^  ^  ^  - 


+  2  k^B^^(0+) 


-  R^i0+)  + 


2  2  2 
^ 


]■- 


Rj (0+)  2 


^ 


)(5-2) 


B„(0+) 
m  * 


-»-  2  kS  ^(0+) 

K. 


2  22 

*1  -  '^  K 


=  0  . 


In  Chapter  IV  it  was  shown  that  as  long  as  the  arch  is 
stable  ^  »  where  m  is  an  unexcited  mode.  Assuming 

that  the  same  will  hold  true  in  this  case  and  then  proving 
that  the  assumption  is  correct  allows  equations  (5*2)  to  be 
treated  as  a  special  case  of  equations  (5*1)  and  so  only  one 
analysis  is  required.  Thus  setting  ~  ^  equations 

(5.1)  gives 


-  -  A/  +  S(T^)j  =  -  R^(t^)  +  k^  -  S^(t^) 


-  .  A/  +  S(x„ 


3) 


-  k^  -  r\2  +  s(t„)]  = 


0 


m  /  l,r 


As  before,  these  equations  should  only  be  used  to  deter¬ 
mine  the  unbuckled  equilibrium  values  of 

B  (t  )  if  they  exist,  under  the  assumption  that  )  =  0. 

r  n  ^  m  n 

The  system  can  then  be  checked  at  t  =  to  determine 

v/hether  B  (t  )  =  0  is  a  valid  assumption  under  the  loading 
m  n 

Rl(-'n)  ' 

Ine  f irr two  of  equations  (5.3)  then  become 


- 

/  +  S(t„)]  =  -  R,(t„)  f  -  S,(t„) 


2  2 
-  k  -  rk 
1 


(5.4a) 


(1  +  E*(T^,T„,^))(r2  +  +  rV(^n)) 


-  +  S(Tjj)  =  -  +  r^(A^  -  S^(Tj^))  .  (S-itb) 

By  maintaining  the  terms  common  to  the  two  equations  on  the 
the  left,  equations  (5.4a)  and  (5.4b)  yield 


T (t  )  +  At  -  S, (t  )  /  . 

i — s - 1 - i 2_  -  (l  +  E  (t  .T  ,  )) 

B,(X)  ^  '  n'  n-lV 

1^  n' 


.  2./. 

+  r  ;^A^ 


Thus  B„(t  )  can  be  expressed  in  terms  of  R,  (t_)  , 
r '  n  I  n 
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-  R,  (t)  +  A,  -  S,  (t„)  ,  ,  ,,  s 

B,(t^)  -l)(l  -^E’t-n'^n-l))  • 

(5-5) 

Substituting  this  result  back  into  equation  (5.^a)  to  remove 
Br(T^)  gives 


J  ^  -  SrKOK(^n) 

+  I  - E - i - 

V-  ^  -1)(1  + 

-  -  v^k^  +  S(Tj^)  =  -  Ri(Tj,)  +  -  S^i-c^)  ,  (5.6) 

a  fifth-order  polynomial  of  • 

Equation  (5.6)  is  plotted  in  Fig.  ^.1  as  a  function  of 
Rl(Tn)  >  wnich  assumes  some  functional  dependence  of 


on  L  J  unbuckled 

equilibrium  value  exists  under  the  assumptions  made  so  far, 
and  the  stable  B^(t^)  vs.  relation  would  be 

given  by  the  curve  ab  of  Pig.  5.I,  which  represents  the 


largest  root  of  equation  (5.6). 


B  (t  )  can  then  be  deter- 
r'  n' 


mined  by  substituting  this  value  of  B^(t^)  into  equation 

(5*  5) • 
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S(V)  =S(t„) 

Si(r„+)  =  S,(t„)  + 

S^(t„+)  =  S^(tJ  +  E*(t„,t„_^)B^(t„) 
Sm('rn'^)  =0  m  /  l,r 


>  (5-7) 


Substituting  these  results  into  equations  (3.12),  as  was 
done  before,  and  by  noting  that  at  t  =  it  is  desired  to 
apply  an  additional  load  AR^  and  AR^  ,  it  follows  that 


B 


=  -  (r^(tJ  +  -  S,(t^+) 

(nir  )  +  AR  ) 

\  r^  n^  r/  , 


r 


r‘-4  -  S  (t  +)) 

\r  r'  n  '/ 


>(5.8) 


■j- 


+  S(T„)J  =  0 


m  /  1,  r 


It  is  convenient  to  consider  AR  as  a  function  of  AR. 

r  1 

so  that  stability  may  be  discussed  with  AR^  the  lone  inde¬ 
pendent  loading  parameter.  While  this  may  not  always  be 


d( 

pi 


tl 

f: 

tv 

I 


an 


Eq 

T7-5 
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desirable,  it  seems  a  reasonable  assumption  to  make  at  tnis 
point . 

In  the  investigation  of  the  system  of  equations  (3.8) 
there  are  again  two  possible  solutions  to  consider.  The 
first  will  be  to  assume  =  0  )  which  allows  the  first 

two  of  equations  (5-8)  to  be  rewritten  so: 


+  B,^(t  +)  + 
1  '  n 


j(. 


A  -  S  (t  (t  +) 

r  r^  n  y  l'^  n  ^ 


-  (Rj^(t„)  +  -  S^(t/)  +  (r2  -  l)B^(T„y 

=  -  (%(^n)  +  -  Si(V)  , 


-  ^  s(t/) 


and 


+  r 


B  (t  ■^)  = 
r '  n  ' 


(5-9) 


—  . (5-10) 


B,  (t  +) 
1 '  n 


+  {r  -  1) 


Equations  (5.9)  and  (5. 10)  are  plotted  as  the  full  lines  on 

T7*?  rro  C  Oo  ^ v%/^  CT 

*  -*.0^  •  <UII\A  ^0C^W  i  ooi-rcv  l/X  • 

The  broken  line  cr  in  Fig.  5.2a  is  given  by 


-  ^  ^1)  ^  ft,  -  Si(t„) 


1  -  r 


(5.11) 
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which  is  found  by  setting  the  denominator  of  equation  (5*10) 

equal  to  zero.  It  can  be  seen  from  equation  (5'9)  that  no 

solutions  of  will  fall  on  this  line  except  at  the 

singular  point  c  ,  and  thus  all  possible  stable  equilibrium 

values  fall  below  this  line  on  the  curve  ab  .  Thus  it  is 

concluded  that  R,  (t  )  +  AR,  =  R^(t  +)  is  the  non-transi- 

tional  buckling  load.  [Rn(T  )]  as  shown  on  Fig.  ^.1  and 

ji.  ri 

Rg(Tn+)  on  Pig.  5.2a  are  equal  only  if  =  f^i^'*^n^^cr  ’ 

If  the  system  is  to  have  an  unbuckled  equilibrium  value  at 

T  =  T  then  it  is  necessary  that 
n 


The  other  possible  solution  to  the  system  of  equations 
(5.8)  has  one  value  of  0  ,  m  f  1,  r  .  Thus  the 
third  of  equations  (5*8)  gives 


+  S(Tj^+)  =  0  . 


(5-12) 

With  the  use  of  equation  (5.12)  the  first  two  of  equations 


(  5.8)  become 


,  (5.13) 
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Substituting  these  two  back  into  equation  (5.12)  results  in 


2„  2 

in 


n 


(5-15) 


If  r  =  2  ,  m  >  2  and  so  the  line  given  by  equation 
(5.13)  would  not  intersect  the  curved  segment  ab  of  Fig. 
5.2a,  since  its  slope  is  greater  than  that  of  the  line  cr  . 
Thus  it  can  be  concluded  that,  if  r  =  2,  ^  ^ 

not  a  possible  solution  and  the  mode  of  buckling  will  be  non- 
transitional  (i.e.  no  unexcited  mode  will  be  present  in  the 
collapse  configuration). 

If  I*  >  3  j  then  m  may  take  on  a  value  less  than  r 

and  it  may  be  possible  to  get  a  solution  where  B  (t  +)  /  0 

m  n 

is  valid.  Equations  (5*13)^  (5*1^)  and  (5-15)  plotted 
as  the  dashed  lines  on  Pigs.  5.2a,  5.2b,  and  5.2c  respective¬ 
ly,  on  the  oasis  of  r  3  and  m  <  r 

As  v/as  the  case  for  the  sinusoidal  arch  with  sinusoidal 
load,  if  b‘  falls  on  the  arc  ab  then  the  critical  loading 
is  ^  inspection  of  equation  (5.13)  it  is  seen 

that  m  =  2  would  give  the  lowest  value  of  R  (t^"^)  .  Un- 
like  the  sinusoiaal  case  though,  it  is  not  convenient  to 
give  Rg(T^'*')  and  in  explicit  form,  since  this 


2 


-62- 


involves  the  solution  of  a  fourth  order  polynomial  rather 
than  a  second. 

The  conclusions  to  be  drawn  are: 

1.  If  Rn  (t  )  <  R„(t^+)  or  R,,(t  ■^)  ^  whichever  is 

X  n  s  ri  u  n 

smaller,  then  the  arch  is  stable  at  ®l^^n^ 

and  found  from  equations  (5.4a)  and  (5.4b),  v/ith 

B  (t  )  =  0  ,  are  the  stable  equilibrium  values, 
m'  n' 

2.  Since  B  (t  )  =  0  as  long  as  the  arch  is  stable, 

m  n 

Sm(Tn)  =  0  ,  and  the  assumption  made  at  the  beginning  of  the 
chapter  that  this  is  so  is  then  proven  to  be  valid. 

3.  If  the  second  mode  is  excited,  then  any  higher  mode 
not  excited  remains  zero  for  all  time,  stable  or  unstable. 

It  is  interesting  to  note  the  influence  of  an  excited 
nonsinusoidal  mode  upon  the  buckling  load.  If  the  second 
mode  is  the  one  that  is  excited,  then  a  very  small  amount 
of  asymmetry  in  the  arch  considerably  lowers  the  critical 
load.  Pig.  5.3  shows  a  typical  load  deflection  relation  at 
tv/o  different  times,  the  dotted  curve  represents  a  sinusoidal 
arch  while  the  full  lines  correspond  to  the  same  arch  except 
that  the  second  mode  has  been  excited.  It  can  be  seen  that 
there  is  considerable  difference  in  the  critical  load  for 
the  two  cases,  and  that  the  influence  of  this  nonsinusoidal 
excitation  increases  with  time. 

Pig.  5.4  shows  the  response  of  the  same  basic  arch  as 
in  Pig.  5.3  except  that  this  case  has  the  third  mode  excited 
instead  of  the  second. 
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Prom  Pigs.  5*3  and  5.4  it  can  be  seen  that  the  effect 
of  exciting  the  second  mode  on  the  buckling  load  of  the  arch 
is  much  greater  than  exciting  the  third  mode.  It  can  also 
be  shown  that  exciting  higher  modes,  either  odd  or  even,  does 
not  have  as  much  effect  as  exciting  the  second  mode. 

As  the  magnitude  of  the  nonsinusoidal  excitation  decreases, 
the  load  deflection  curves  approach  the  dotted  curve  of  the 
sinusoidal  arch  and  the  cr  line.  Thus  in  the  limit  as  the 
excitation  goes  to  zero,  the  critical  load  for  the  nonsinu¬ 
soidal  arch  approaches  the  critical  load  for  the  sinusoidal 
arch,  as  would  be  expected. 

Pigs.  5.3  and  5.4  correspond  to  arches  that  have  an 

initial  rise-to-span  ratio  such  that  the  mode  of  buckling  is 

2  /  +  \ 

transitional,  which  requires  that  >5*5  ^ 

where  t  is  the  time  at  which  the  arch  buckles.  If 
cr 

p  , 

A^  -  <  5.5  ^  the  sinusoidal  arch  will  buckle  in  a 

nontransitional  manner.  Pig.  5.5  shows  the  response  of 
such  an  arch  with  higher  modes  excited. 

It  can  be  seen  from  Pig.  5^5  that,  for  A^  =  A^  ,  the 
effect  of  the  second  mode  being  excited  is  again  larger  than 
the  effect  due  to  the  third  mode. 


Comparison  of  Pigs.  5.3  and  5.4  with  5.5  show  that  for 


constant 


the  influence  of  the  nonsinusoidal  mode 


on  the  buckling  load  decreases  as  A^  decreases.  Por  the 
elastic  case  [11 ]  shows  plots  of  critical  load  vs.  A^  for 
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various  ratios  of 


and 


and  these  shov;  the  above 


mentioned  trend  very  clearly. 


> 


(V.II)  Several  modes  excited 

The  solution  of  the  problem  is  which  many  modes  are  ex¬ 
cited,  through  either  the  initial  arch  ahape  or  the  loading, 
is  a  simple  extension  of  the  case  of  two  modes  excited. 

Equations  (5.I)  through  (5*5)  remain  valid  except  that 
r  now  stands  for  any  excited  mode  except  the  first  mode. 

Equations  (5.6)  can  then  be  written  as 


+  2  k 


(1  + 


r 


5 


rC“-r 


/ 

=  -  ,  (5-16) 


a  polynomial  of  of  order  (2h  +  1)  ,  v^here  h  is  the 

total  number  of  excited  modes. 

Equation  (5.16)  is  plotted  in  Pig.  5.6.  The  dashed  curve 
in  Pig.  5.6  is  given  by 


(1 


2  2 
-  2  ir  A. 

k=r  ^ 


+  S(r„) 


- 


+  A. 


-(£- 
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6 


<7 


and  it  can  be  seen  by  comparing  this  to  equation  (5.16)  that 
all  real  solutions  of  equation  (5.16)  must  fall  within  the 
dashed  curve  as  shovm. 


If  Rt (t  )  <  (t  )  then  an  unbuckled  equilibrium 

J-  “  Li  n  j  cr 

position  exists  and  is  represented  by  a  point  on  the  arc  ab 
Bi(Tn)  can  be  calculated  from  equation  (5.16)  and  the  other 
modes  then  calculated  from  equation  (5* 5)-  Pollov/ing  the 
same  procedure  as  was  used  when  only  two  modes  were  excited 
S(Tn+)  ,  and  =  0  can  be  calcu¬ 

lated,  and  equations  (5.8)  then  give  the  governing  system  of 
equations  for  the  time  t  =  . 

Again  two  possible  solutions  exist,  the  first  with 

B  (t  +)  =0  which  gives 
m  n 


B.  (t  +) 
I  n 


1  +  B  ^(t  ■^) 
1  '  n  ^ 


+  2  k 
k=r 


)  +  ^J  0/  \\ 

.(RlUn> 


-  +  S(t„+) 


=  -  +  S,(tJ  , 

(5.17) 


and  as  given  by  equation  (5. 10). 

The  second  solution  with  one  ^ 

Bi(t^)  and  B^(t^+)  as  given  by  equations  (5-13)  and  (5.14) 
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respectively,  and  governed  by 


where  m  ^  l,r 

If  the  second  mode  Is  excited,  then  as  before  B  (t  +) 
will  always  remain  zero,  but  if  the  second  mode  is  not  ex¬ 
cited,  ^  ^  where  m  =  2  may  be  a  possible  solution. 

Fig.  5*7  shows  equations  (5*13)  and  (5.17)  for  the  case 
m  =  2,  i.e.  the  second  mode  is  not  excited.  Thus  if  b' 
falls  along  the  portion  ab  as  shown,  criti¬ 

cal  loading,  but  if  b'  falls  along  be  or  does  not  exist 
then  R  (t  +)  is  critical. 

It  can  be  seen  that  the  only  difference  in  the  solution 
of  the  case  where  many  modes  are  excited,  as  compared  to  the 
case  of  two  excited  modes,  is  that  the  order  of  the  poly¬ 
nomial  that  determines  B^(t)  ,  equations  or  (5.17)^ 

is  raised  by  two  for  every  additional  excited  mode.  While 
this  raises  the  possibility  of  having  more  equilibrium  posi¬ 
tions,  it  is  noted  that  these  additional  positions  are  not 
unbuckled  stable  equilibrium  positions,  and  so  in  the  actual 
evaluation  of  the  critical  buckling  load  they  are  not  a  factor. 
The  fact  that  the  polynomial  is  of  higher  order  causes  some 
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increase  in  the  work  required  to  calculate  the  unbuckled 
equilibrium  position,  but  this  is  the  essential  difference 
between  the  two  cases. 

Pigs.  5‘8  and  5.9  show  the  ordinate  B^(t)  plotted 
against  r  ,  the  dimensionless  time  parameter,  for  a  fluid 
and  a  solid  material  respectively.  Aside  from  the  fundamental 
difference  in  the  two  figures,  caused  by  the  fact  that  for 
loads  below  a  certain  value  the  arch  made  of  a  solid  material 
reaches  a  stable  equilibrium  value  as  t  ,  the  response 

of  the  two  cases  is  qualitatively  the  same.  It  is  to  be 

noted  that  for  small  values  of  asymmetry  <  0.001 

^1 

the  ordinate  is  nearly  identical  with  that  of  the  sinusoidal 
case  until  the  critical  condition  is  approached,  at  v/hich 
time  it  starts  to  diverge  rather  rapidly  with  subsequent 
failure . 

Pigs.  5.10  and  5. 11  are  plots  of  the  critical  buckling 
time  vs.  the  applied  load  for  various  amount  of  asjrmmetry  in 
the  second  mode.  Of  particular  interest  is  the  large  decrease 
in  the  critical  time  for  very  small  amounts  of  this  asymmetry. 
As  pointed  out  previously  the  effect  of  exciting  the  third 
or  higher  mode  is  not  as  great  as  exciting  the  second  mode, 
but  the  character  of  the  response  is  the  same  as  that  shown 
in  Pigs.  5.8  to  5.11. 
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i.me 


Three  Parameter  Solid  Material 
Sinusoidal  Load  At  =  4 


Chapter  VI  Long-time  Solution  for  a  Viscoelastic  Solid 
Material 

If  the  arch  is  made  of  a  viscoelastic  solid  material, 
and  the  loading  approaches  a  constant  as  t  becomes  large, 
then  there  is  a  possibility  of  the  deflection  reaching  some 
final  stable  value  as  t  . 

Letting  Tp  denote  t  as  it  becomes  very  large,  equa¬ 
tions  (3.21)  can  be  written  as 


.  2f.  /7PE'(Tp-i) 


m"  ‘P 


pP  E  /  /  \  \  1 

'K  'L  -Ef— 


(6.1) 


m  =  1,2,3, ...  . 

In  the  evaluation  of  the  integrals  it  is  to  be  noted  that 
B  (i)  is  essentially  constant  and  equal  to  B  (t„)  for 
^  Tp  .  For  i  small  varies  but  E‘(Tp-4)  0 

since  (xp-^)  becomes  very  large.  Therefore  in  the  range 
of  i  where  E'(t^-4)  is  different  from  zero 
nearly  constant  and  so  can  be  taken  out  from  under  the  in¬ 


tegral  sign.  Thus  to  a  good  approximation 


-’■p  E'(T„-e) 


E„  -  E, 


Where  Ep  =  E(<»)  . 


(6.2) 
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Similarly, 


(6-3) 


Substituting  equations  (6.2)  and  (6.5)  into  equations 
(6.1)  gives 


o  “f"  m  A 


,(6.4) 


m  =  1,2,  ... 

Equations  (6.4)  are  identical  with  the  equations  of  an 

elastic  ai^h  of  modulus  under  the  same  loading;  but 

in  the  viscoelastic  case  they  can  only  be  used  to  determine 

the  unbuckled  equilibrium  position  denoted  by  the  ordinates 

B  (t„)  ,  since  in  their  derivation  it  was  assumed  the  arch 
m  F 

was  stable  under  the  loading  R  (t„)  and  had  reached  a 

rn  r 

final  configuration. 

To  check  the  stability  at  t  =  Xp  ,  subject  the  arch 
to  an  additional  load  at  x  =  Xp"^  in  the  same  manner  as 
was  done  previously  in  Chapters  IV  and  V.  B^(xp)  can  be 
calculated  from  equations  (6.4),  considering  any  unexcited 
modes  to  remain  zero.  If  an  unbuckled  solution  of  equations 
(6.4)  does  not  exist,  this  indicates  the  arch  has  buckled 
at  some  previous  finite  time.  On  the  assumption  that 
Bj^(xp)  exists,  the  integrals  S(xp'*')  and  S^(xp‘'')  (see 
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The  solution  for  this  set  of  equations  and  the  method 

of  determining  stability  are  given  in  Chapter  IV  or  Chapter 

V,  depending  on  the  number  of  excited  modes .  If  it  is  found 

that  (AR. )  <  0  ,  this  indicates  that  the  arch  has  buckled 

at  some  previous  finite  time. 

It  is  to  be  noted  that  the  critical  load  found  from 

equations  (6.6)  differs  from  the  critical  load  of  a  purely 

elastic  arch  of  modulus  that  has  the  same  deflected 

o 

shape  as  the  long-time  solution  of  the  viscoelastic  arch 
under  the  same  loading  •  'fhis  is  due  to  the  fact 
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that,  although  the  two  arches  support  the  same  load  v/ith  the 
same  deflected  shape,  the  axial  thrust  and  moment  distribu¬ 
tion  are  not  the  same  in  both  cases.  This  arises  because  the 
moment  is  directly  proportional  to  the  deflection  while  the 
thrust  is  proportional  to  the  square  of  the  displacement. 

In  [11 ]  it  was  shovm  that,  for  a  uniform  load  on  a  sinu¬ 
soidal  arch,  the  effect  of  the  higher  modes,  i.e.  3^  5^  !> 
v;as  to  reduce  the  critical  load  by  no  more  than  .  5^  of  the 
critical  sinusoidal  load.  That  this  is  not  valid  for  the 
viscoelastic  case  can  be  expected  from  the  results  of  Chapter 
V,  where  it  is  seen  that  near  the  critical  buckling  time  the 
displacement  of  the  non-sinusoidal  arch  differed  considerably 
from  the  displacement  of  the  sinusoidal  arch  (see  Fig.  5.9) • 
The  reason  for  this  can  be  seen  from  an  examination  of  equa¬ 
tions  (5.17)  and  (5.18).  In  each  of  these  expressions  there 
is  a  series  of  the  form 


v/here  r  includes  any  excited  mode  other  than  the  first. 

In  order  to  illustrate  the  effect  of  the  viscoelastic  materi¬ 
al,  consider  the  case  of  a  uniform  load  on  a  sinusoidal  arch 

which  has  R-,  (t) 

Rp('t’)  =  )  r  =  3j  5^  7 f  ••• 

=0  ,  r=2,4,6,...  , 

and  -  0  ,  r  /  1  . 
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Substituting  these  into  the  series  (6.7)  results  in 


2 

.(6.8) 


k= 


,  A  2  2\2 

=3^ 5w  •  •  •  -  m  j 


In  the  elastic  case  S^(t) 
duced  to 


k=-  c; 


.7. 


0  and  the  series  (6.8)  is  re- 


(6.9) 


v.’hlch  is  a  very  rapidly  converging  series.  For  m  =  2 


2 

k=3,5w 


k 


-  .y 


=  4.977  X  10 


-4 


and  so  the  effect  of  the  higher  modes  Oi.  the  elastic  critical 
load  is  seen  to  be  negligible.  Kov/ever,  for  the  viscoelastic 
case  ^  ^  ^  second  term  of  the  series  (6.8) 

leads  to  an  expression  of  the  form 


2 

k=3,5.7 


(6.10) 


which  could  be  divergent,  depending  on  the  values  of  S  (t  ■*■)  . 

r  n 

While  it  is  not  possible  to  express  explicitly  S.^(t^'^)  as 
a  function  of  r  ,  it  is  seen  that  in  the  long-time  solution 
S^(Tj^‘^)  varies  directly  as  B^(Tp)  (see  equation  (6.5)  ), 
and  that 
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4 


(6.il) 


E,,  Hj^(Tp) 


Br(Tp)  = 


E. 


_ Q 


+  (r  -  1) 


r  =  3.  5,  7 


Equations  (6.11)  are  derived  from  equations  (6.4)  for  a 
uniform  load  on  a  sinusoidal  arch.  Thus  it  is  seen  that 


for  large  r  ,  varies  as  — ^  and  so  the  conver¬ 

gence  of  the  series  (6.10)  is  assured.  The  fact  remains 
however,  that  the  series  (6.8)  may  not  be  negligible  since 
some  of  the  lower  terms  may  have  appreciable  magnitude,  and 
cannot  be  neglected  in  the  calculations.  Added  to  this  is 
the  fact  that  the  terms  S(Tp+)  and  Sj,(Tp+)  appear  in 
other  ways  in  equations  (5.17)  and  (5.18)  and  must  be  ac¬ 
counted  for. 

Pigs.  6.1  to  6.4  show  the  effect  of  Ep  and  R^(Tp) 
on  the  buckling  load  of  an  arch  with  initial  rise  =  4  . 
Compared  to  Fig.  6.1,  which  is  the  purely  sinusoidal  case, 
the  other  figures  show  a  very  rapid  decrease  in  the  critical 
load  as  R2(Tp)  increases.  This  can  be  attributed  to  the 
fact  that  the  amplitude  of  the  higher  modes  in  the  unbuckled 
configuration  increases  very  rapidly  as  the  load  approaches 
the  critical  load  (see  Pig.  5*2). 
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cr 


Ep  =  -7 


.5 


pT.  12 


Fig.  6.1 


T’h'ree  Parsipstsr*  Solid  j  t  t 

Sirusoidal  Arch  A,  =  4  ,  Sinusoidal  Load 


Chapter  VII  General  Features  of  the  Numerical  Solution 


The  numerical  integration  and  solution  of  the  roots  of 
the  algebraic  equations  were  performed  on  the  B-5OOO  digital 
computer  at  Stanford  University. 

In  the  numerical  integration  of  integrals  of  the  type 

T 

f(e)de  (7.1) 

it  is  advantageous  to  use  a  variable  time  step  if  the  form 
of  the  functions  E(t)  and  f(T)  are  of  the  requisite  type. 
E’(t)  varies  quite  rapidly  for  t  small  but  as  t  becomes 
large  E’(t)  approaches  zero.  In  this  problem,  where  fCr) 
represents  functions  of  the  type  ”  "^m  )  ^ 

possible  to  draw  several  conclusions  regarding  the  behavior 
of  f(T)  if  the  loading  functions  sufficiently 

smooth.  For  t  small  f(T)  will  vary  quite  rapidly,  re¬ 
gardless  of  the  type  of  viscoelastic  material,  but  for  t 
large  the  material  property  makes  a  considerable  difference. 
For  a  fluid  material  it  is  expected  that  f(T)  would  con¬ 
tinue  to  vary  quite  rapidly  until  failure  takes  place,  es¬ 
pecially  if  the  time  at  failure  is  quite  small.  On  the  other 
hand  if  the  material  is  a  viscoelastic  solid  the  variation 
at  large  times  will  be  quite  small,  except  near  the  critical 
time . 

In  order  to  take  advantage  of  these  facts  and  reduce 
the  computing  time,  a  logarithmic  time  step  increment  was 


used  (see  [I5])  until  the  critical  time  was  found,  at  which 
time  the  program  reverted  back  to  the  last  known  stable  time 
and  a  smaller  time  step  introduced.  In  this  manner  small  time 
steps  were  used  for  t  small  and  t  near  the  critical  time, 
while  larger  steps  were  used  where  E‘(t)  and  fCr)  varied 
more  slov/ly. 

Among  the  roots  >  only  the  largest  root  of  B^(t) 

is  required  in  the  analysis  and  so  a  simple  iterative  type 
solution  was  used.  Convergence  was  very  rapid  for  solutions 
near  the  point  a  of  Pig.  4.1  but  became  slovier  as  the  equi¬ 
librium  value  approached  b  ,  at  vjhich  point  the  program  would 
not  converge . 

If  the  progrmm  would  not  converge  a  check  was  made  to  in¬ 
sure  that  no  unbuckled  equilibrium  solutions  existed,  as  this 
constitutes  one  possible  mode  of  failure.  The  second  possible 
type  of  failure,  that  of  an  unexcited  mode  becoming  non-zero, 
is  easily  checked  by  calculating  on  the  basis  of 

/  0  (see  for  example  equation  (4.9)  )/  then  if  the 
B^(t)  found  in  this  manner  is  greater  than  the  B^(t)  found 
from  the  iterative  solution,  the  arch  is  unstable. 

The  time  step  increments  and  the  allowable  tolerance  in 
the  iterative  solution  were  adjusted  so  that  the  solution  for 
B^(t)  was  accurate  to  at  least  four  significant  digits  in 
all  the  test  cases  considered.  The  maximum  computer  time 
needed  was  less  than  one  minute  for  any  one  problem  including 
a  case  where  ten  modes  were  excited  and  the  maximum  time  taken 
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to  be  1000. 


In  Chapter  VI  there  are  no  time  integrations  to  be  per¬ 
formed,  but  in  order  to  find  tlie  critical  load  for  equations 
(6.6)  it  was  necessary  to  use  a  double  iterative  procedure. 
One  iterative  procedure  is  required  to  determine  the  un¬ 
buckled  equilibrium  position,  and  a  second  is  then  required 
to  find  the  critical  buckling  load.  Thus,  although  the 
programming  is  quite  simple,  considerable  machine  time  is 
required  in  its  execution. 
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